First-principles study of the incorporation and diffusion of helium in cubic zirconia 
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The incorporation and difTusion of helium (He) with and without intrinsic vacancy defects in cu- 
bic Zr02 are investigated through first-principles total-energy calculations, in which the projector- 
augmented-wave (PAW) method with the generalized gradient approximation (GGA) is used. The 
calculated formation energies of intrinsic point defects indicate that cubic Zr02 has a tolerant resis- 
tance to radiation damage. The incorporation energy of He impurity shows that it is preferable to 
occupy the Zr vacancy at first, whereas the solution energy suggests that He would be accommodated 
in the interstitial site at thermodynamic equilibrium concentration. By calculating the He migra- 
tion energies corresponding to both interstitial and vacancy assisted mechanisms, we suggest that 
it is most likely for He to diffuse by hopping through a single vacancy. Remarkably, our calculated 
vacancy-assisted diffusion energy of He is consistent well with the experimental measurement. 
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Zirconia (Zr02) is a highly attractive material of great 
interests both scientifically and in terms of its techno- 
logical applications due to its excellent mechanical, ther- 
mal, chemical, and dielectric properties For example, 
Zr02 is employed as superplastic structural ceramic that 
demonstrates superb strength and fracture toughness , 
and is also used as one of the best refractive engineer- 
ing materials for thermal barrier coating on aeronauti- 
cal and land-based gas turbine blades Q. Additionally, 
the cation-doped Zr02 can also be used in solid oxide 
fuel cells and oxygen sensors Furthermore, as one 
of the most radiation resistant ceramic materials, Zr02 
is applied as a promising candidate host phase for the 
inert matrix nuclear fuel (IMF) in order to immobilize 
and burn up 6] the plutonium from dismantled nu- 
clear weapons and other high level nuclear waste (HLW) 
caused by the actinides (e.g. ^*^Am). Consequently, 
great attentions are needed to focus on the basic scientific 
research of Zr02. 

Zr02 exists at least five crystalline structures with dif- 
ferent symmetries. The monoclinic P2\/c polymorph is 
the only one found at ambient condition whereas the 
tetragonal P^i/nmc and the cubic FmUm phases are sta- 
ble above 1400 and 2600 K and the other two orthorhom- 
bic Phca and Pnma phases are stable above 3 and 20 
GPa, respectively |8l-[l2|. As these five Zr02 structures 
are centrosymmetric, they are all nonpolar and nonpiezo- 
electric. Presently, Zr02 stabilized in the cubic phase 
with e.g. Y2O3 is considered as a suitable refractory 
host for actinide confinement used for transmutation or 
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high-temperature reactor projects. Large amounts of ex- 
perimental measurements have been carried out to de- 
termine the stability of this material in high radiation 
environment. One early investigation on phase stabil- 
ity of Zr02 under fission fragment damage was reported 
by Wittels and Sherrill [l^ . Then the evolution of radia- 
tion induced damages in Zr02 Q was systematically in- 
vestigated using Rutherford backscattering spectrometry 
and ion channeling (RBS/C), along with X-ray diffrac- 
tion and transmission electron microscopy (TEM) un- 
der Xe and I ion irradiation. No amorphization of the 
pure Zr02 was observed. But Meldrum et al. [isj found 
that nanocrystalline Zr02 could be amorphized under 
Ne and Xe ion irradiation. Damage evolution, surface 
morphology and characteristic structural changes in He- 
implanted cubic Zr02 were extensively investigated in 
experiments p^| - [l9| . Fluence dependence and thermal 
stability of defects in He- implanted cubic Zr02 was deter- 
mined by slow positron implantation spectroscopy [2(|. 
Moreover, the He damage and He effusion in fully stabi- 
lized Zr02 [m was measured with neutron depth profil- 
ing (NDP) for different annealing temperatures and ana- 
lyzed by electron microscopy during the various anneal- 
ing stages concluding that the diffusion of He is probably 
caused by vacancy assisted interstitial diffusion. Besides, 
He migration in monoclinic and cubic yttria-stabilized 
Zr02 was also studied from non-destructive '^He depth 
profiling using the resonant '^He(d,p)*He nuclear reaction 

From theoretical point of view, ah initio methods have 
been employed to study the structural, electronic, me- 
chanical, thermodynamical, and vibrational properties of 
Zr02 [U, [23l - [33| . Very recently, the energetics of intrinsic 
point defects and the influence of Al doping in Zr02 [s^ - 
|36| have been studied to identify the dominant defects un- 
der different oxygen chemical potentials and Fermi levels. 
However, up to date there are no detailed investigation 
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performed on the He incorporation and diffusion behav- 
ior in Zr02 through a state-of-the-art first-principles ap- 
proach. To make a deep understanding of the effect of 
He in Zr02 in atomic scale is helpful for us to be more 
confident to validate the application of Z1O2 in severe 
environment. Therefore, this is the main purpose of the 
present work. In the present study, first the behavior of 
He in cubic Zr02 is investigated by determining the for- 
mation energy of point defects. Then the incorporation 
energy of He into different sites of the lattice is calcu- 
lated, through which the most favorable incorporation 
site is determined. By taking thermodynamic equilib- 
rium into account, the solubility of He is also evaluated. 
After exploration of the stability of He, we subsequently 
concentrate on the diffusion mechanism of He in Zr02 by 
calculating energy barriers and migration pathways be- 
tween two trapping sites. After comparison of calculated 
results for different postulated mechanisms, we present 
the viable diffusion mechanism of He in Zr02 matrix. 
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FIG. 1: Conventional unit cell of Zr02. O atoms are located 
at the corners of small cubes, and Zr atoms are located at 
the center of an alternative cubes, x indicates the octahedral 
interstitial site (OIS) in the lattice. 



III. RESULTS AND DISCUSSIONS 



A. Crystallography and bulk properties 



II. COMPUTATIONAL METHOD 



The density functional theory (DFT) total energy cal- 
culations are self-consistently performed in the frame- 
work of the frozen core all-electron PAW method [37] 
as implemented in the Vienna ab initio simulation pack- 
age (VASP) • The electron exchange and correlation 
potentials are evaluated within the GGA introduced by 
Perdew, Burke, and Ernzerhof (PBE) [s^. Twelve elec- 
trons (4s^4p^4d^5s^) for zirconium (Zr) and six electrons 
(2s^2p^) for oxygen (O) are taken into account as valence 
electrons. To get accurate results, the cutoff energy of the 
plane wave expansion is set up to a high value of 800 eV, 
which was tested to be well converged with respect to 
the total energy. The integration over the Brillouin Zone 
(BZ) is performed with a grid of special k point-mesh de- 
termined according to the Monkhorst-Pack scheme (40| . 
After convergence test, 8x8x8 k point-mesh is chosen to 
calculate the bulk properties with the total energy differ- 
ence less than 1 meV per unit cell. In the present work, 
large supercells containing up to 96 atoms have been em- 
ployed to reduce any artificial error due to the use of a 
smaller supercell as it will be shown later. Correspond- 
ingly, a 2x2x2 A;-point grid is used to sample the BZ of 
the 96-atom supercell for the modeling of the point de- 
fects and He impurity in Zr02. Relaxation procedures 
at ground state are carried out according to the Gaus- 
sian broadening with a smearing width of 0.1 eV. For all 
defect structures, the atomic relaxation is considered to 
be completed when the Hellmann-Feynman forces on all 
atoms are less than 0.01 eV/A. And the volume relax- 
ation of the supercell is also calculated. 



Here, we considered the cubic phase of Zr02, which 
is crystallized in the fluorite structure with space group 
Fm3m (No. 225). The Zr atoms are in the center of a 
cube of O atoms, while the O atoms occupy the tetrahe- 
dral sites coordinated by Zr atoms, constituting a sim- 
ple cubic sublattice as shown in Fig. 1. The present 
optimized equilibrium lattice parameter (a) obtained by 
fitting the energy-volume data in the third-order Birch- 
Murnaghan equation of states (EOS) is 5.152 A and 
some other basic bulk properties calculated with this 
PAW method is also listed in Table 1, in fair agreement 
with the experimental and previous theoretical results. 



B. Intrinsic point defects formation 

The previous study of Zr02 bulk properties gives in- 
sights into the reliability of the DFT-GGA approach to 
model the cubic Zr02. Now we turn to the study of the 
stability of different types of point defects in the fluorite 
lattice. The first-principles method allows us to consider 
Zr02 supercells containing up to hundred atoms and thus 
to carefully check the convergence of the defect formation 
energies as a function of the supercell size. Indeed, in 
too small a supercell and due to the periodic boundary 
conditions, a point defect interacts with its image in the 
adjacent supercell, which yields a spurious contribution 
to the calculated formation energy. To check this conver- 
gence, we have considered supercells representing repeti- 
tion of the elementary cubic cell according to the follow- 
ing repetition pattern: 1x1x1 (12 atoms), 2x1x1 (24 
atoms), 2x2x1 (48 atoms), 2x2x2 (96 atoms), 3x2x2 
(144 atoms), 3x3x2 (216 atoms), 3x3x3 (324 atoms). 
The defects are explicitly investigated in the calculation 
of formation energy as follows: vacancies, intcrstitials at 
OIS, Frenkel pair and Schottky defect. 
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TABLE I: Calculated equilibrium lattice parameter oo (A), bulk modulus Bq (GPa), and cohesive energy -Ecoh (eV) of fluorite 
Zr02 at ground state. For comparison, available experimental values and other theoretical results are also listed. 

Parameter Present work Previous calculation Experiment 

ao (A) 5.152 5.130", 5.128", 5.116" 5.086", 5.108% 5.090-' 

Bo (GPa) 234 267", 251% 235"= 194", 215*" 

£coh (eV) -15.99 -23.28% -11.45-' -11.45^ -14.72' 

"Ref. [27], "Ref. [30l,=Ref. [33], "Ref. [42], 
"Ref. [43], ^Ref. [44], "Ref. [45], Ref. [46], 
' Ref. [26], ' Ref. [47], ''Ref. [4^], ' Ref. [49] 



The formation energies and Efp of a vacancy ( Vx) 
and an interstitial (/x) of the X specie are obtained from 
the total energies of the system with and without the 
defect, according to 

E^^ ^ E^^'^- ~ Ezro, + Ex, (1) 

and 

E'f = ^Z?02 ~ ^Zr02 - ^x, (2) 

where i?zi02 the energy of the defect-free Zr02 super- 
cell, E^^Q^ and E/^q^ are the energies of the supercell 
containing one vacancy and interstitial, respectively. Ex 
is the energy per atom of each chemical species in its ref- 
erence state {X= Zr or O). Here the reference state of Zr 
is chosen as the ground state crystalline phases, namely, 
the a-Zr crystal with space group Pd^/mmc (No. 194). 
As for the O, the atomic oxygen energy is employed as a 
reference state [s^] • In the total energy calculation of an 
O atom, an orthorhombic supercell with lattice constants 
exceeding 15 A is employed and the spin-polarization is 
included. The calculated formation energy of the O va- 
cancy in cubic Zr02 is in good agreement with previous 
theoretical and experimental results [13] ■ 

Frenkel pair consists of a non-interacting vacancy and 
an interstitial of the same chemical element. Thus the 
formation energy can be determined from the formation 
energies of the vacancy and of the interstitial calculated 
separately, given by 

i?;^- = + (3) 

A Schottky defect is a more complex defect consisting of 
a Zr vacancy and two O vacancies, all of which are again 
supposed to be non-interacting. The formation energy is 
calculated by 

E'f'' = ElTo. + '2EZ-^^ ~ 3^£;z,.02 , (4) 

where N is the total atom number in the used supercell. 

The variation in the calculated formation energies of 
the single point defects as a function of the size of the 
supercell is represented in Fig. 2. On the whole, clearly, 
one can see that the convergence becomes satisfied when 
the supercell size is chosen as large as 2x2x2 (96 atoms). 
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FIG. 2: Formation energies of single point defects in Zr02 as 
a function of the size of the supercell used in the calculation. 



TABLE II: Calculated formation energies Ef (eV) of point 
defects in Zr02: Zr and O vacancies (Vac), octahedral in- 
terstitials (Int), Frenkel pair (FP), and Schottky defect (Sh). 
The energies are calculated with atom relaxation (pos) and 
with atom and volume relaxation (pos-|-vol). Aa/a is the 
variation in the lattice parameter of the 96-atom supercell. 



Parameters 


Int Zr Int O Vac Zr Vac O FP Zr FP 


Sh 


Ef (eV) pos 


2.14 1.19 15.27 8.60 17.40 9.79 


18.68 


Ef (eV) pos-l-vol 


1.97 1.17 15.26 8.60 17.22 9.77 


18.66 


Aa/ a 


0.5% 0.1% 0.1% -0.1% - 


0.2% 



Therefore, all the results reported in the rest of the paper 
are obtained with the 2x2x2 supercells. Table 2 gives 
the formation energies of Zr and O vacancies, interstitials 
at OIS, Frenkel pair, and Schottky defect, all of which are 
obtained after performing the atom relaxation, and with 
or without volume relaxation. The variation in the lattice 
parameter Aa/a of the supercell is also reported. This 
variation is only indicative of the degree of convergence 
of the size of the supercell: ideally, for an infinite-size 
supercell, this variation should be zero. 

It can be seen from Table 2 that the formation ener- 
gies of various point defects are positively large, consis- 
tent with the high resistance property of Zr02 against 
radiation damage. Furthermore, Table 2 shows that the 
O defects have much lower formation energies than Zr 
defects. Thus radiation damage and deviation from sto- 
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ichiometry will be preferably accommodated in the O 
sublattice: by the formation of O vacancies in hypostoi- 
chiometric Zr02 and by the formation of O interstitials 
in hyperstoichiometric Zr02, for instance. The effect of 
volume relaxation on the formation energies is negligi- 
ble, except for the Zr interstitial for which it is 0.2 eV. 
For this defect, the supercell size variation is the largest, 
but does not exceed 1%. Unfortunately, to our knowl- 
edge, the experimental values of the formation energies 
of point defects for cubic Zr02 are not available in the 
literature for comparison. 



C. Stability of He impurity in Zr02 



TABLE III: Incorporation Einc (eV) and solution energies i?soi 
(eV) of He in Zr02 at different sites of the fluorite lattice: 
the OIS, the Zr and the O substitution sites with relaxation 
of the atom position (pos) and volume of the supercell (vol), 
change in the supercell lattice parameter Aa/a, and the elec- 
tron charge density Q of He at the trap sites. 



Parameters 



He(OIS) He(Zr) He(0) 



-Bine (eV) pos 
-Bine (eV) pos+vol 
-Bsoi (eV) pos 
-Bsoi (eV) pos+vol 
Aa/a (%) 
Q (electrons) 



1.3491 
1.3457 
1.3491 
1.3457 
0.03 



0.4332 1.5931 
0.4199 1.5779 
7.1327 6.4871 
7.1176 6.4650 
0.15 0.12 



2.0596 2.0298 2.1891 



The stability of He at pre-existing trap sites is de- 
termined by calculating the incorporation energy which 
is defined as the energy needed to trap a He at a pre- 
existing trap site (vacancy or interstitial considered here) 
in the expression as follows: 



and 



-Bine — 



He 
Zr02 



(5) 



where E^°q_^ is the energy of Zr02 supercell containing 
the He impurity, Ezr02 is the energy of the supercell 
with an empty trap site of He, and -Ehc is the energy of 
an isolated He. Therefore, a positive value means that 
the energy is required to incorporate He in the lattice 
whereas a negative result implies that incorporation is 
energetically favorable. A comparison of incorporation 
energies is the most simple way to assess the stability of 
the impurity and predict the most stable solution site for 
He provided that trap sites are available for occupation 
[Slj . This will be the situation when the concentration 
of He atoms is low enough that incorporation proceeds 
through occupation of intrinsic defect sites. 

The use of incorporation energy is limited since it is 
not sensitive to any equilibrium between trap sites. To 
take into account the equilibrium between the different 
trap sites one should consider solution energies. Solu- 
tion energies should be used to express the population of 
He in the different sites when complete thermodynam- 
ical equilibrium is achieved. For a certain site X, the 
solution energy of He is defined as the incorporation en- 
ergy plus the apparent formation energy of the trap site. 
It is clear that for an interstitial site incorporation and 
solution energies are e qual . For vacancy insertion sites, 
the expression is given [S^l as 



Eu 



E 



(6) 



va- 



where E^^ is the apparent formation energy of the 
cancy site. With stoichiometry at the ground state, the 
apparent formation energies of Zr and O vacancy sites 
can be expressed 52] as 



E 



^FPo 



(7) 



2 ^ 



(8) 



Here Ep'' and EfJ^ are the O Frenkel pair and Schottky 
defect formation energies. 

We investigate the He behavior at various locations in- 
cluding the high-symmetry positions in the lattice, and 
find that He always moves to the OIS during atomic re- 
laxation in a defect-free Zr02. For defective structures, 
we take considerations of trap sites with one Zr or O 
vacancy as the substitution site, respectively. The incor- 
poration energies and solution energies are calculated by 
taking into account the relaxation of the atomic positions 
in the supercell, with and without the volume relaxation, 
and for the different incorporation sites. In addition, 
through Bader analysis [H, HI] , we have also calculated 
the electron charge density of He when it occupies var- 
ious trap sites. All the relevant results are presented in 
Table 3. 

The calculated incorporation energies at any trap site 
are positive indicating that they are less stable in Zr02 
fluorite lattice than in the free gas phase. This means 
that some energy has to be provided to incorporate He 
in Zr02. The minimum incorporation energy is found 
for the Zr vacancy. With relatively high incorporation 
energies of He at OIS and O vacancy, we expect a slight 
preference for He to be located at Zr vacancy. However, 
the large solution energies at various sites (see Table 3) 
lead to the prediction of the insolubility of He in cubic 
Zr02. We can conclude that the driving force for the 
experimentally observed He release has its origin in this 
insolubility. The results shows that He is more likely 
to be accommodated at the OIS for the thermodynamic 
equilibrium concentration. As for the perturbation of 
the crystal induced by the impurity, it is revealed that 
the volume relaxation with the incorporated He has no 
great influence on the incorporation energies and the He 
impurity only induces a slight expansion of the lattice pa- 
rameter. Finally, Bader analysis results show that the He 
impurity displays only a weak charge transfer, in agree- 
ment with a complete electronic shell and neutral envi- 
ronment at this position. This finding also explains why 
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a Zr vacancy is the most favorable incorporation site for 
He since it is well known [s^ that a He atom prefers to 
occupy the site with the low electron density due to its 
fiUed-shell electronic configuration. 

D. Diffusion of He impurity in Zr02 

The concentration of He remaining in the lattice is de- 
cided not only by the solution energies but also by the 
ease with which the He atom can diffuse out of the lat- 
tice. Here, we take considerations into the possible diffu- 
sion mechanisms for one He atom. The climbing image 
nudged elastic band (CINEB) method [s^ is employed 
to find the minimum energy path (MEP) of atomic mi- 
gration when both the initial and final configurations are 
known. The CINEB is an efficient and improved method 
for finding the energy saddle point between the given 
initial and final states of a transition to understand the 
energy path of He diffusion behavior. Migration energy 
barrier is defined as the energy difference between an 
initial configuration and the saddle point. The intersti- 
tial and vacancy assisted diffusion mechanisms of one He 
atom in Zr02 are systematically studied as follows. 

1. interstitial diffusion mechanism 

We first investigate the interstitial diffusion mechanism 
of He in defect-free Zr02 by calculating the migration en- 
ergy between two adjacent OIS sites, which are the body- 
centered position of the lattice. The diffusion pathway 
and MEP are shown in Fig. 3. We perform the CINEB 
calculation in which the initial guess for the atomic posi- 
tions of the diffusing He atom is simply a linear interpo- 
lation between the OIS sites. The saddle point obtained 
with this initial guess is found to be in the middle of the 
migration pathway as shown in Fig. 3(b), corresponding 
to an energy barrier of 3.74 eV. In the panel of Fig. 3(b), 
we can see obviously the distortion of the nearest two O 
atoms at the saddle point, which are pushed away from 
their lattice sites by about 0.37 A in the perpendicular 
direction to the migration pathway of He by the strain 
energy. 

2. Zr vacancy assisted mechanism 

Next, we explore the migration pathways of He in a 
defective Zr02 and find that the energy barriers are re- 
markably decreased compared to the previous results. 
We take Zr and O vacancy into account separately to 
probe into the minor difference of the vacancy assisted 
diffusion mechanism. In the case of Zr vacancy preexist- 
ing, three possible migration pathways are assumed and 
CINEB calculations are carried out similarly. 

The first postulated mechanism is involved in the di- 
rect diffusion of He from OIS into the nearest OIS in the 
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FIG. 3: (a) Diffusion pathway and (b) MEP of He migration 
from OIS into the nearest OIS in defectfree Zr02. 



presence of Zr vacancy, shown in Fig. 4(a). Examining 
the MEP for this mechanism in Fig. 4(b), two images 
on either side of the saddle point configuration exhibit a 
slight oscillation in energies around the transition state. 
The amplitude of the oscillation is only 0.02 eV, which is 
within the error of the calculations. The reliable estimate 
of the energy barrier can be obtained as about 0.23 eV, 
which must be overcome in the direct diffusion path of He 
from OIS to OIS. We find that the Zr vacancy is a trap 
for He atom from the MEP. Along the diffusion pathway. 
He atom moves into the Zr vacancy during the atomic 
relaxation at the third image, and an additional energy 
of about 1.36 eV must be provided when the He migrates 
back into the nearest OIS as shown in Fig. 4(b). The dis- 
tortion of the nearest O atoms of He is calculated. Due 
to the strain energy of He, the displacement of O atom 
in the perpendicular direction to the pathway is about 
0.43 A at the saddle point. 

The second kind of diffusion is assumed that He hops 
through a Zr vacancy from one OIS to the next nearest 
OIS. In addition, the two adjacent OIS lies in the line 
with Zr vacancy, represented in Fig. 5(a). It is interest- 
ingly noted that the MEP in Fig. 5(b) displays the same 
migration energy trend as the first one in Fig. 4(b). The 
energy barrier is 0.23 eV for migration from OIS to the 
nearest OIS, and once the He atom diffuses into the Zr 
vacancy, it hardly moves into the next OIS unless an ex- 
tra energy of 1.36 eV is gained. The next nearest O atom 
has the largest distortion 0.48 A when He is trapped in 
the Zr vacancy. 

The third possible diffusion path is similar with the 
second assumption, but the two nearest OIS is not in 
the line with Zr vacancy, as illustrated in Fig. 6(a). We 
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FIG. 4: (a) Diffusion pathway and (b) MEP of He migration 
from OIS into the nearest OIS in the presence of Zr vacancy 
in Zr02. 
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FIG. 6: (a) Diffusion pathway, (b) MEP of He migration from 
OIS into Zr vacancy, and (c) MEP of He migration from the 
Zr vacancy into the nearest OIS. 



tion energy barrier than the interstitial diffusion mecha- 
nism. Different diffusion pathways proposed here demon- 
strate a similar migration characteristic. Obviously, it is 
more likely to trap the He atom into the Zr vacancy along 
the diffusion pathway when compared to the defect-free 
cases. 



FIG. 5: (a) Diffusion pathway and (b) MEP of He migration 
from OIS into the next nearest OIS in the presence of Zr 
vacancy in Zr02. 



consider the CINEB calculation of this migration in two 
steps. The MEP in Fig. 6(b) shows the step that the He 
atom first jumps into the Zr vacancy from the OIS, and 
Fig. 6(c) gives the situation that He atom migrates into 
the nearest OIS from Zr vacancy. From initial OIS to 
Zr vacancy, the saddle point is obtained with an energy 
barrier of 0.23 eV, and when He diffuses from Zr vacancy 
to final OIS, another saddle point appears with a larger 
migration barrier of 0.91 eV. Similarly, in this case the 
distortion of the nearest O atoms of He at Zr vacancy is 
as small as only 0.10 A. 

From above analysis, it is concluded that in Zr02, as Zr 
vacancy preexists. He diffuses easily with a lower migra- 



3. O vacancy assisted mechanism 

As for O vacancy that preexists in Zr02, similar con- 
sideration is taken to evaluate the diffusion mechanism 
of He atom. At first, we discuss that He atom directly 
migrates from OIS into the nearest OIS, given in Fig. 
7(a). Then He atom hopping from OIS to nearest OIS 
through O vacancy is investigated. It is considered that 
the He atom first moves into the O vacancy from the 
OIS, and next the He atom is assumed to jump into the 
nearest OIS from the O vacancy, shown in Fig. 8(a). The 
MEP of direct migration from OIS to the nearest OIS in 
the presence of O vacancy is different from the migra- 
tion with the preexisting Zr vacancy discussed above, as 
shown in Fig. 7(b). There are two images in Fig. 7(b) 
that are predicted to be most energetically favorable. In 
the middle point of the pathway, the metastable state is 
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FIG. 7: (a) DifTusion pathway and (b) MEP of He migration 
from OIS into the nearest OIS in the presence of O vacancy 
in Zr02. 



found between the two local saddle points. At the middle 
image of this pathway, the He atom is trapped into the O 
vacancy after the relaxation during CINEB calculation. 
The energy barrier from the energy favorable site to the 
saddle point is 2.18 eV. On the other hand, the MEP of 
He atom hopping from OIS to nearest OIS through O 
vacancy is illustrated separately in Fig. 8(b) and Fig. 
8(c). The saddle point is obtained in the pathway with 
the energy barrier 1.04 eV and 0.92 eV, corresponding 
the migration energy barrier from the initial OIS to O 
vacancy and from O vacancy to the final OIS, respec- 
tively. At various images in the CINEB calculation, the 
O atoms near the O vacancy move inward to the vacancy 
at about 0.25A. 

From the above discussion, we find that the migra- 
tion barrier of interstitial diffusion is relatively higher 
than the vacancy assisted diffusion. During the migra- 
tion of He with Zr or O vacancy preexisted, the He atom 
is easily trapped into the vacancy and an extra energy 
of only about 1 eV is needed to hop out of the vacancy, 
which remarkably agrees well with the experimental re- 
sult [531 . Therefore, the diffusion behavior of He atom in 
cubic Zr02 is dominated by the vacancy assisted diffusion 
mechanisms. 



IV. CONCLUSION 




I Mr ■ ^ 

if- 1 1^ 



» Zr • O ;> VacO • He X OIS 
(b) 




FIG. 8: (a) Diffusion pathway, (b) MEP of He migration from 
OIS into O vacancy, and (c) MEP of He migration from the 
O vacancy into the nearest OIS. 



that trapping of a He atom at a Zr vacancy is more sta- 
ble than at a O vacancy or at an OIS. From the solution 
energy results, on the other hand, it has been found that 
the He is more likely to be located at the OIS for the 
thermodynamic equilibrium concentration. To get more 
insights into the diffusion mechanism of He, we have cal- 
culated the migration energy of He in various pathways 
through the CINEB calculation, from which it has been 
shown once the He is trapped by the vacancy, its mi- 
gration energy barrier from the vacancy to the nearest 
OIS is largely decreased when compared to the vacancy- 
free case. Thus, we conclude the diffusion of He in Zr02 
is mainly assisted by the vacancy, which is expected to 
provide a guiding line in explaining the experimentally 
observed He diffusion phenomenon. 



In summary, we have performed total-energy calcula- 
tions to investigate the formation energies of native point 
defects, the He incorporation and solution energies, and 
the diffusion properties of He in cubic Zr02, using the 
PAW-GGA method and supercell approach. Our cal- 
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